Inhomogeneity and the Post-Inflationary 

Universe 



Richard Easther and Matthew Parry 



Q\ . Department of Physics, Brown University, Providence, RI02912 

Abstract. We discuss the interaction between perturbations in the inflaton and the 
background during the preheating phase in simple inflationary models. By numerically 
00 ■ solving the Einstein field equations we are able to assess the impact of non-linear grav- 

es) , itational effects on preheating, and to assess the accuracy of perturbative discussions 

of the preheating epoch. 

> 

cn ■ INTRODUCTION 

^ Inflation drives the primordial universe towards the very special initial state 

needed for it evolve into the form it is observed to have today. However, to be suc- 
cessful, inflation must terminate gracefully with the energy density of the inflaton 
field being transferred to radiation, a process known as reheating. 

Since 1990 [1] it has been realized that reheating can be driven by nonlinear, 
coherent effects leading to explosive particle production via parametric resonance. 
The resulting distribution is far from thermal equilibrium, and the process is often 
called preheating (see [2,3], and refs. within). 

The standard analytic approach to preheating is based on a Fourier expansion 
of the relevant fields. The equation of motion of the k-th mode has the form of 
a Mathieu or Lame equation. Modes corresponding to values of k in resonance 
bands grow exponentially, and the precise resonance structure depends sensitively 
on the underlying particle physics. Because it selectively amplifies specific Fourier 
modes, parametric resonance implies an enhancement in the spatial variation of the 
fields. While the non-linear backreaction on the field evolution has been examined 
analytically [2,4] and numerically [5,6] the backreaction of the inhomogeneity on the 
underlying spacetime metric and its implications for the evolution of the universe 
has only recently begun to be studied [3,7-13]. 

We have broken new ground by tackling the back-reaction problem using the 
full Einstein field equations, and eschewing the use of any perturbative approxima- 
tions. We show that large metric inhomogeneities can be induced by parametric 
resonance, so approximations to the full field equations may not adequately de- 
scribe the evolution of the universe during (and after) preheating. At present, our 



principal simplifying assumption is that the inhomogeneity lies in a single spa- 
tial direction, which allows us to work with a 1 + 1 dimensional system of partial 
differential equations. 

Here we briefly survey the growth of metric perturbations during preheating in 
an inflationary model driven by a X(p 4 potential. When treated analytically, the 
m 2 4> 2 model does not undergo parametric resonance during reheating [11], a result 
we confirmed numerically [3] . However, A0 4 possesses a single resonance band, and 
when we examine the evolution of modes within this band we see observe strong 
amplification of the corresponding metric perturbations. This work will be more 
fully described in our forthcoming paper [13]. 



METRIC AND INITIAL CONDITIONS 

We have assumed that the universe has a planar symmetry, so the metric func- 
tions depend only on t and z, and are independent of x and y. In [3], we studied 
reheating after m 2 (j) 2 inflation using the metric 

ds 2 = dt 2 - A 2 (t, z) dz 2 - B 2 (t, z) (dx 2 + dy 2 ), (1) 

which describes an inhomogeneous universe in which the dx 2 + dy 2 sections have 
zero spatial curvature. In principle, we could retain this metric for the A</> 4 case, 
but it is advantageous to work in conformal-like co-ordinates analogous to those 
which simplify the homogeneous system [4,13]. Specifically, by transforming t and 
z to r] and ( we write equation (1) as: 

ds 2 = a 2 (7], C) (dr] 2 - d( 2 ) - (3 2 (t], () (dx 2 + dy 2 ). (2) 

In practice we use the co-ordinate transformation to compute A and B from equa- 
tion (1), as well as a and (3, as we evolve the system numerically. The equa- 
tions of motion will be given in [13], and the numerical techniques used to solve 
them are similar to those we used in [3]. We focus on initial configurations where 
a single mode (when viewed in the perturbative context) is excited. We choose 
a(0,C) = 13(0,0 = 1, 0(0,0 = 0o and a c (0, () = fl c (0, C) = 0,f(O,C) = 0. The 
constraints are solved if 




2 " } + V(<m) = C. (3) 



«,,(o, = ^ + vw) - §, A„(o, = ^1 

Our choice of C ensures (a tV (0X)) = (Ar?(0, C)), where (• • •) is a spatial average. 
The actual inhomogeneity is injected through the inflaton kinetic energy, 

^(0,C) = H + esin(^), (4) 

where Z is the length of our "box", and n is the average initial velocity. Since 
both t and z are transformed, slices of constant rj in the conformal frame are not 
mapped directly to slices of constant t in the physical frame. The initial density 
perturbation is on the order of e 2 . 
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FIGURE 1. The evolution of A (the g zz component of equation (1)) is plotted as a function of 
r\ and £. The left panel shows the evolution of a mode with k slightly too large to be in resonance, 
while the right hand case undergoes resonance, and significant inhomogeneity is generated. The 
units are arbitrary. The initial perturbation is e 2 = 10~ 6 , and the simulations begins at the end 
of inflation. 

RESULTS 

When we examined reheating after m 2 4> 2 we confirmed that there is no significant 
resonant amplification, and that the perturbative analysis is valid [3] . However, af- 
ter X(p A inflation there is a single, narrow resonance band, and the scaling properties 
of the solution ensure that modes which are initially in this band remain there in- 
definitely. 

Fig. 1 shows the evolution of the metric component A for two different pertur- 
bations, one inside the resonance band and one outside it. A significant degree of 
inhomogeneity is generated in the metric by the resonance in the fields. While a 
simple analytical treatment suggests that resonance lasts indefinitely, the backreac- 




FIGURE 2. The Fourier mode of $, the metric perturbation, for a resonant mode is shown: 
the left panel gives the perturbative result, while the right panel shows the evolution of the mode 
derived from the full nonlinear analysis. 



tion from the nonlinear field evolution eventually halts the resonant amplification. 
This feature is apparent in our simulation, and is illustrated in Fig. 2, where we plot 
the metric perturbation, $ [14]. The termination of resonance coincides roughly 
with the comparatively sudden growth of the metric perturbation seen qualitatively 
in the right panel of Fig. 1. 

If we solve for an inhomogeneous field, <fi, while averaging over the spacetime 
background (corresponding to numerical treatments which include the nonlinear 
field equations, but assume that the metric is unperturbed) the subsequent evo- 
lution of $ differs from that which we obtain in the presence of the gravitational 
backreaction. Comparing the two results will allow us isolate the contribution from 
the inhomogeneous parts of the metric to the overall evolution. We are also ex- 
tending this work to more general initial perturbations than equation (4), so we 
can include interactions between different modes and study non-perturbative effects 
such as the formation of primordial black holes after inflation [13]. 
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